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Abstract 

This  article  presents  an  adaptive  backstepping  approach 
to  the  longitudinal  (7,  a,  Q)  control  of  an  aircraft  that 
directly  accommodates  magnitude,  rate,  and  bandwidth 
constraints  on  the  aircraft  states  and  the  actuator  sig¬ 
nals.  The  article  includes  design  of  the  control  law,  Lya¬ 
punov  analysis  of  the  stability  properties  of  the  closed 
loop  system,  and  simulation  based  analysis  of  the  per¬ 
formance. 

1  Introduction 

The  introduction  of  uninhabited  combat  air  vehicles 
(UCAV’s)  has  generated  interest  in  adaptive  control  for 
aircraft  that  are  capable  of  maintaining  trajectory  fol¬ 
lowing  control  even  subsequent  to  faults  or  battle  dam¬ 
age.  In  addition,  physical  limitations  impose  magni¬ 
tude,  rate,  and  bandwidth  constraints  on  the  control 
surface  deflections  and  the  aircraft  states.  In  this  article, 
we  present  analysis  and  simulation  results  for  an  adap¬ 
tive  backstepping  approach  that  is  designed  to  accom¬ 
modate  magnitude,  rate,  and  bandwidth  constraints  on 
the  actuator  signals  as  well  as  the  intermediate  control 
variables  used  in  the  backstepping  control  approach. 

Adaptive  control  has  been  applied  to  aircraft  in,  for  ex¬ 
ample  [2,  3,  4,  14,  15,  18].  Piloted  hardware  in  the  loop 
simulation  and  flight  testing  is  described  in  [2,  4,  18]. 
The  approach  of  [3,  4]  used  a  dynamic  inversion  based 
nonlinear  controller  with  a  neural  network  approximat¬ 
ing  the  error  in  the  inversion  process.  Adaptive  back- 
stepping  has  been  applied  to  aircraft  in  for  example  [15]. 
That  article  did  not  consider  the  effects  of  saturation  on 
the  closed  loop  performance. 

Control  signal  rate  and  amplitude  constraints  in  adap¬ 
tive  control  are  addressed  in,  for  example,  [1,  7,  9,  10, 
17] .  The  first  type  of  method  is  to  completely  stop  adap¬ 
tation  under  saturation  conditions.  This  ad-hoc  method 
does  prevent  the  tracking  errors  induced  by  actuator 
constraints  from  corrupting  the  parameter  adaptation, 
but  is  undesirable  as  it  does  not  allow  any  theoreti¬ 
cal  stability  guarantees.  The  second  set  of  approaches 
(training  signal  hedging  (TSH)),  see  e.g.  [1,  9],  cor¬ 
rects  the  tracking  error  only  where  it  is  used  in  the  pa¬ 
rameter  adaptation  laws.  It  does  not  directly  change 


the  reference  input  to  the  control  loop.  The  third  set 
of  of  techniques,  referred  to  as  pseudo-control  hedging 
(PCH)  [7,  10],  only  alters  the  commanded  input  to  the 
loop.  The  goal  in  PCH  is  to  decrease  the  command  to 
the  loop  to  a  point  where  it  is  implementable  without 
saturation.  PCH  does  not  directly  change  the  signal 
used  in  the  parameter  adaptation  laws.  The  work  in  [7] 
develops  and  analyzes  the  PCH  method  to  extend  the 
control  approach  of  [3,  14]  to  accommodate  magnitude 
and  rate  limits  on  the  commanded  surface  deflections. 
Although  the  second  and  third  approaches  start  with 
radically  different  philosophies,  for  first  order  systems, 
it  can  be  shown  that  these  two  approaches  are  identical. 
For  higher  order  or  nonlinear  systems,  the  comparison 
is  not  as  straightforward. 

This  article  presents  an  adaptive  backstepping  approach 
for  longitudinal  control  of  an  aircraft  that  has  state  and 
actuator  constraints.  Both  theoretical  and  simulation 
analyses  are  included.  A  novel  aspect  of  this  approach  is 
the  ability  to  accommodate  magnitude,  rate,  and  band¬ 
width  constraints  on  the  actuator  signals  and  each  of  the 
intermediate  control  variables  (i.e.,  aircraft  state  vari¬ 
ables)  used  in  the  backstepping  design.  This  result  is  of 
interest  in  its  own  right.  This  approach  has  similarities 
with  both  the  TSH  and  PCH  methods,  since  it  corrects 
both  the  training  signal  and  the  reference  input  at  each 
step  of  the  backstepping  procedure.  Our  ultimate  goal 
is  to  extend  this  approach  to  a  full  vehicle  controller 
using  on-line  function  approximation  [5,  11,  12]  to  ac¬ 
commodate  battle  damage  types  of  events. 

2  Problem  Formulation 

Due  to  space  limitations,  the  discussion  of  this  article 
will  focus  on  a  simplified  aircraft  model  with  the  same 
structure  as  the  linearized,  rigid  body,  aircraft  longitu¬ 
dinal  dynamics  [16].  The  model  is  just  complex  enough 
to  allow  a  presentation  and  analysis  of  the  control  de¬ 
sign  approach.  Since  the  model  has  the  same  structure 
as  aircraft  longitudinal  dynamics,  the  approach  can  be 
extended  to  a  realistic  nonlinear  aircraft  model.  This 
extension  is  discussed  further  in  the  conclusions  section. 
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The  model  that  we  will  consider  has  the  form 

7  =  La  +  Laa 
ol  =  Q  —  ( La  +  Laa) 

Q  =  M0  +  MqQ  +  M5S, 

where  the  (normalized1)  lift  force  is  modeled  as  L0+Laa 
and  the  (normalized2)  pitch  moment  is  modeled  as 
Ma  +  MqQ  +  Mg 5.  In  this  expression,  7  is  the  flight 
path  angle,  a  is  the  angle  of  attack,  Q  is  the  pitch  rate, 
and  5  is  the  control  surface  deflection.  At  an  operating 
point,  La,  La,  Ma,  Mq  and  Ms  are  unknown,  constant, 
scalar  parameters.  Additional  terms  can  be  added  to 
the  lift  or  moment  models  without  changing  the  theo¬ 
retical  approach.  Most  importantly,  the  approach  di¬ 
rectly  extends  to  the  case  where  multiple  independent 
surfaces  are  available.  These  generalizations  are  not  in¬ 
cluded  here  due  to  space  limitations. 

The  corresponding  tracking  error  dynamic  equations  are 


3  Adaptive  Backstepping  with  Saturation 

The  backstepping  control  laws  are  summarized  below. 
The  operating  envelope  requires  a  to  remain  in  the  in¬ 
terval  [«£,  a^]  and  Q  to  remain  in  the  interval  [Q^ ,  Q(?]. 
The  actuator  design  requires  5  to  remain  in  the  interval 
[5c,  5%\.  In  addition,  the  actuator  design  will  impose 
rate  and  bandwidth  limitations  on  <5,  which  in  turn  will 
result  in  rate  and  bandwidth  limitations  on  the  inter¬ 
mediate  variables  Qc  and  ac.  The  above  operating  en¬ 
velope  and  surface  deflection  constraints  are  assumed  to 
be  known  and  fixed. 

The  nominal  backstepping  commands  are 

(~^°  +  'Vc  -  7  +  ?y7)  -  xa 

Qc  =  ( Lo  +  Lao^j  -  kaa  +  +  ac  -  7 La  -  \Q 

5c  =  -  (~M0  —  MqQ  —  fcqQ  +  r/Q  —  a  +  QCJ 


7  =  La  +  La  (ac  +  a)  -  7 c 
a,  =  —  (La  +  Laa)  +  Qc  +  Q  —  ac 

Q  —  Ma  +  MqQ  +  Ms5  —  Qc 


with  fc7,  ka,  kQ  >  0.  The  robustifying  terms  (??7,  r]a,  ijq) 
are  motivated  in  the  proof  of  the  theorem  to  follow.  The 
amplitude,  rate,  and  bandwidth  limited  control  signals 
can  be  defined  by 


where  the  tracking  errors  are  defined  as  7  =  7  ^  7C, 

xc>  Q  —  Q  Q 


ea  =  (sat(a0c,otc,aLc)  -  xai) 


a  =  a  —  OLr 
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will  cause  7  to  track  an  input  signal  7C.  We  assume  that 
7 c  is  smooth  and  that  7 c  is  known.  In  the  backstepping 
approach,  the  7  controller  will  generate  a  commanded 
value  for  a.  An  a  controller  will  generate  a  Q°c  com¬ 
mand  for  Q.  The  <5  controller  will  generate  the  control 
surface  command  5°  to  achieve  the  Q  command. 
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Since  an  aircraft  operating  envelope  includes  constraints 
on  both  actuators  and  states,  the  controller  must  work 
in  the  presence  of  magnitude,  rate,  and  bandwidth  lim¬ 
its  on  5  as  well  as  operating  envelope  constraints  on 
Q  and  a.  In  the  formulation  to  follow,  we  will  treat 
(L0,La,M0,MQ,Ms)  as  parameters  to  be  identified. 
Therefore,  we  are  implementing  adaptive  backstepping 
with  control  variable  and  state  constraints.  In  the  ac¬ 
tual  application,  this  set  of  linearization  ‘parameters’  is 
a  function  of  the  operating  conditions.  Since  the  oper¬ 
ating  condition  is  slowly  time  varying,  the  adaptive  laws 
would  be  required  to  adapt  the  parameters  to  match  the 
local  model.  In  the  final  full  vehicle  implementation, 
these  ‘linearized  parameters’  would  be  approximated  as 
a  function  of  the  flight  condition. 


1  To  simplify  the  notation,  the  factor  at  the  operating  point 
is  included  in  the  lift  function. 

2To  simplify  the  notation,  the  -M  factor  is  included  in  the 

1yy 

moment  function. 


where  the  parameters  of  this  filter  are  defined  below.  In 
this  filter,  ea  is  the  error  between  the  first  filter  state 
and  the  magnitude  limited  filter  input.  Since  the  first 
state  of  the  filter  is  used  in  the  computation  of  ea,  the 
characteristic  equation  of  the  filter  is  s2  +  2C,aLjnas  + 
io\a  =  0,  in  the  linear  portion  of  the  two  sat  functions. 
The  sat  function  in  the  second  equation  enforces  the 
rate  limit.  At  the  filter  output,  ac  is  magnitude,  rate, 
and  bandwidth  limited.  The  command  filters  for  Qc 
and  5C  are  defined  similarly  as 


eQ  =  (sat  (Q°cQc,Qc)  -  xQl) 
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e5  =  (sat  (8°, 6^,6%)  -  xSl) 
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In  these  expressions,  La  =  2<f%Ra ,  Lq  =  2’%Rq  , 

na  nQ 

Lg  =  2<^Rs  and  (Ra,  Rq,  Rg)  are  rate  limits  on  (a,  Q ,  <5), 

n<5 

respectively;  (LOna,LOnQ,uns)  are  bandwidth  limits  on 
0 a,Q,S );  (a^,Q^,S^)  are  upper  limits  on  (a,Q,8); 
(af,QR,5f)  are  lower  limits  on  (a,  Q,  S),  and  sat  de¬ 
notes  the  saturation  function  that  is  linear  with  unity 
slope  between  its  lower  and  upper  limits.  This  para¬ 
graph  has  presented  one  command  filter  approach, 
many  other  filters  are  possible  without  affecting  the  va¬ 
lidity  of  the  stability  proofs.  If  the  surface  position  can 
be  measured,  then  the  last  filter  is  not  required.  Also 
the  parameters  of  the  above  filters  are  not  completely 
independent.  For  example,  since  the  derivative  of  a  is 
closely  related  to  Q ,  the  rate  limit  Ra  should  be  less 
than  or  equal  to  the  magnitude  limit  Q(f . 

The  vector  x  =  [x7(i), Xa(*), XqW]  is  generated  by 

X7  -  ^7X7  T  ^OL  («C  ^c) 

Xa  ^aXtt  “b  ( Qc  Qc) 

XQ  =  -kQXQ  +  Mg  (8C  -  5°)  . 


These  signals  are  filtered  versions  of  the  effect  of  state 
and  control  rate,  bandwidth,  and  magnitude  constraints 
on  the  variable  that  is  being  controlled.  Finally,  we  de¬ 
fine  the  vector  of  modified  tracking  errors  x  =  [7,  a,  Q\ 
as 


control  rate,  bandwidth,  and  magnitude  constraints.  In 
an  implementation,  deadzones  should  be  used  to  avoid 
parameter  drift  due  to  noise;  and,  projection  methods 
must  be  used  to  ensure  that  the  estimation  transient 
does  not  cause  La  or  Mg  to  change  sign. 

4  Tracking  Error  Dynamics 

Representing  the  actual  parameters  as  the  estimated  pa¬ 
rameter  minus  parameter  estimation  error:  La  =  La  — 
L0,  La  =  La  -  La,  M0  =  M0  -  tiQ,  Mq  =  Mq  -  Mq, 
Mg  =  Mg  —  Mg;  the  dynamics  of  the  modified  tracking 
errors  are 


7 

=  7  -  x7 

=  ( La  +  Laot 

-7c)  - 

^  ^7X7 

+  Lq 
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-Mo- 

mqQ- 

MgS 
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(6) 

The  final  numbered  equations  will  be  used  in  the  sub¬ 
sequent  stability  analysis. 

5  Lyapunov  Analysis 

The  objective  of  this  section  is  to  prove  the  stability 
of  the  closed  loop  system  using  Lyapunov-like  methods. 
The  results  are  summarized  in  the  following  theorem. 

Theorem  5.1  The  closed  loop  system  with  dynamics 
described  in  Section  2  and  controller  described  in  Section 
3  has  the  following  properties: 


7  =  7  -  X7 
a  =  a  —  Xa 
Q  =  Q~xq- 

When  there  is  no  control  variable  saturation,  x  con¬ 
verges  to  x.  In  the  presence  of  state  and  control  vari¬ 
able  constraints,  the  signal  x  is  designed  to  remove  the 
effects  of  the  saturation  from  x  to  generated  x,  which 
will  be  used  in  the  parameter  adaptation  laws.  The 


parameter  adaptation  laws 

are: 

i'R  =  r30 

L0  =  (7-  a) 

(1) 

iiQ  =  t4qq 

La  =  r2  (7  -  of)  « 

(2) 

tig  =  r  5qs, 

(3) 

where  T,  for  i  =  1, . . . ,  5  are  positive  adaptation  gains. 
The  subsequent  sections  will  show  that  the  above  con¬ 
trol  laws  are  stable  even  in  the  presence  of  state  and 


1.  L0,La,M0,MQ,Mg  £  £00 

2.  7,  a,  and  Q  £  C <*> 

3.  7 ,  a,  and  Q  converge  to  zero  with  7,  a,  and  Q  £ 
C  2. 

In  addition,  although  we  do  not  show  it  here,  it  is  possi¬ 
ble  to  show  that  given  certain  persistence  of  excitation 
conditions,  the  parameter  errors  converge  exponentially 
to  zero. 

Proof:  Define  the  7-Lyapunov  function  as 

v=-2(f  +  a2  +  Q 2)  +  ^r-1© 

where  0T  =  [L0,  La,  Ma,  Mq,  Mg],  0  is  the  estimate  of 
0,  and  0  =  0  —  0.  Taking  the  time  derivative  of  V, 
using  eqns.  (4-6)  to  eliminate  (7,  a,  Q),  and  using  eqns. 


(1-3)  to  eliminate  0,  after  some  algebraic  manipulations 
yields 

V  =  - kQQ 2  -  kaa2  -  fc7 y2 

+7?y7  +  0:770,  +  Qt]q.  (7) 

Therefore,  the  designer  selects  (rjj,  rja,  t]q)  as  odd,  dif¬ 
ferentiable  functions  with  sign  opposite  to  (7,  a,  Q),  re¬ 
spectively.  The  magnitude  of  these  robustifying  terms 
would  be  selected  to  ensure  boundedness  of  tracking  er¬ 
rors  in  the  presence  of  model  error,  measurement  noise, 
or  disturbances.  For  the  results  of  Section  6,  all  three 
of  these  terms  are  set  to  zero. 

Since  V(t)  is  nonincreasing,  we  conclude  that 
(<3,d,7,  0T)  are  each  bounded.  Barbalat’s  Lemma  im¬ 
plies  that  the  origin  of  the  ( Q ,  a,  7)  variables  is  asymp¬ 
totically  stable.  Finally,  because  0  <  V(t)  <  V(0),  inte¬ 
gration  of  both  sides  of  (7)  shows  that  (Q,  a,  7)  £  C-2-^) 

This  theorem  shows  that  even  when  the  nominal  com¬ 
mands  (a°c,Q°c,50)  exceed  the  state  and  control  limita¬ 
tions,  the  quantities  x  =  [7 ,a,Q]T  and  0  do  not  di¬ 
verge.  In  fact,  the  cumulative  norm  of  these  quantities 
as  defined  by  y/V  ( t )  will  never  increase.  Nothing  is 
proved  about  the  tracking  error  x,  which  may  increase 
during  periods  when  the  state  or  control  limitations  are 
in  effect,  because  the  desired  control  signal  is  not  being 
implemented  (i.e. ,  uc  7^  it°).  The  vector  %  is  defined  to 
remove  from  x  the  portion  of  the  tracking  error  that  is 
due  to  violation  of  the  state  or  control  signal  limitations. 
The  compensated  error  x  is  used  to  adapt  the  parame¬ 
ters,  even  during  periods  in  which  the  state  or  control 
signal  limitations  prevent  implementation  of  the  desired 
control  signals.  When  the  state  or  control  signal  limi¬ 
tations  are  not  in  effect,  (i.e.,  uc  =  tt°),  %  approaches 
zero,  and  x  converges  towards  x. 

6  Simulation  Example 

This  section  analyzes  results  from  a  simulation  example. 
The  parameters  for  the  simulation  were  defined  in  Table 
1.  The  control  gains  are  selected  to  enforce  a  time  scale 
separation  between  the  three  loops  of  the  backstepping 
controller.  The  section  contains  two  sets  of  results.  The 
first  set  of  results  pertain  to  a  situation  where  saturation 
does  not  occur.  The  second  set  of  results  pertain  to  a 
situation  where  significant  saturation  occurs  in  all  three 
control  loops.  The  parameters  of  the  command  filters 
defined  in  Section  3  are  Ca  =  Cq  =  wBil  =  ka,  and 

X>nQ  =  kQ. 

6.1  Example:  No  saturation 

For  this  set  of  simulations,  7r  is  a  25s  square  wave  with 
peak  magnitude  of  ±5  deg.  The  commands  yc  and  jc 
are  the  states  of  a  second  order,  relative  degree  two, 
unity  gain  prefilter  with  (  =  1,  and  u>n  =  fc7. 

Figure  1  contains  plots  of  the  states  and  actuator  signal 
for  the  first  and  last  25  seconds  of  a  150  s  simulation. 


Each  graph  contains  three  curves.  The  dotted  curve  is 
the  command.  The  solid  curve  is  the  magnitude,  rate, 
and  bandwidth  limited  command.  The  dashed  curve 
is  the  actual  state  or  actuator  variable.  Even  though 
the  initial  parameter  error  is  large,  the  transient  in  the 
control  response  is  reasonable.  Figure  2  includes  graphs 
of  the  three  components  of  x  in  the  left  column  and  the 
three  components  of  x  in  the  right  column.  Note  that 
as  the  theory  predicted,  since  there  is  no  saturation, 
x  — >  x  — >  0. 

The  convergence  of  x  to  zero  is  not  monotonic.  How¬ 
ever,  the  Lyapunov  function  is  non-increasing.  The  left 
column  of  Figure  3  plots  the  norm  of  the  parameter  er¬ 
ror,  the  norm  of  x,  and  the  norm  defined  by  the  square 
root  of  the  Lyapunov  function.  Notice  that  although 
the  norm  of  x  is  both  increasing  and  decreasing,  the 
Lyapunov  function  never  increases.  The  right  column 
of  Figure  3  plots  each  component  of  0. 

6.2  Example:  Substantial  Saturation 

For  this  set  of  simulations,  7,.  is  a  25s  square  wave  with 
peak  magnitude  of  ±10  deg.  The  commands  7 c  and  jc 
are  the  states  of  a  second  order,  relative  degree  two, 
unity  gain  prefilter  with  £  =  1,  and  u>n  =  fc7. 

Figure  4  contains  plots  of  the  states  and  actuator  signal 
for  the  first  and  last  25  seconds  of  a  150  s  simulation. 
Even  though  the  initial  parameter  error  is  large,  the 
transient  in  the  control  response  is  reasonable.  Note 
that  the  large  7  commands  result  in  rate,  magnitude, 
or  bandwidth  constraining  the  aci  Qc,  and  8C  variables. 
Figure  5  includes  graphs  of  the  three  components  of  x 
in  the  left  column  and  the  three  components  of  x  in  the 
right  column.  Note  that  as  the  theory  predicted,  since 
there  is  saturation,  x  — >  0,  but  x  does  not  converge  to 
zero.  In  spite  of  this  fact,  the  parameter  estimates  con¬ 
verge  toward  the  correct  values.  In  fact,  as  shown  in 
the  right  column  of  Figure  6,  the  parameters  converge 
faster  than  they  did  in  Figure  3.  There  is  a  significantly 
higher  level  of  excitation  in  the  present  example,  but 
the  parameter  estimates  will  only  converge  if  the  train¬ 
ing  error  is  properly  compensated  to  remove  the  effects 
of  the  rate,  magnitude,  and  bandwidth  limitations  on 
the  control  signals.  The  left  column  of  Figure  6  contains 
plots  of  the  norm  of  the  parameter  error,  the  norm  of 
x,  and  the  norm  defined  by  the  square  root  of  the  Lya¬ 
punov  function.  Note  that  the  Lyapunov  function  is 
non-increasing  even  during  the  periods  of  saturation. 

7  Conclusions 

This  article  has  presented  and  analyzed  an  adaptive 
backstepping  approach  to  the  control  of  longitudinal 
aircraft  dynamics  that  directly  accommodates  magni¬ 
tude,  rate,  and  bandwidth  limits  on  the  aircraft  states 
and  actuators.  For  this  preliminary  analysis,  we  have 
only  considered  a  linearized  aircraft  point  design.  Due 
to  the  positive  analysis  of  the  performance,  we  are  cur- 


rently  extending  this  research  in  four  directions.  First, 
we  are  extending  the  adaptive  state  and  actuator  con¬ 
strained  backstepping  control  approach  to  include  con¬ 
trol  of  lateral-directional  dynamics.  Second,  we  are  ex¬ 
tending  the  linearized  adaptive  approach  to  a  nonlin¬ 
ear  on-line  function  approximation  based  approach  [6]. 
Third,  we  are  extending  the  theoretical  analysis  as  a 
general  extension  of  the  backstepping  approach[13].  Fi¬ 
nally,  we  have  extended  the  results  presented  herein  to 
the  case  where  the  vehicle  has  multiple,  redundant  and 
independent  actuators. 

8  Acknowledgement 

This  research  supported  by  University  of  California- 
Riverside,  and  University  of  Cincinnati,  and  Barron  As¬ 
sociates  Inc.  through  Air  Force  PRD  A  01-03- VAK. 

References 

[1]  A.  M.  Annaswamy  and  J.  E.  Wong.  “Adaptive  control 
in  the  presence  of  saturation  nonlinearity”,  Int.  J.  of  Adap¬ 
tive  Control  and  Signal  Processing,  Vol.  11:  3-19,  1997. 

[2]  J.  Brinker  and  K.  Wise,  “Flight  testing  of  a  recon- 
hgurable  flight  control  law  on  the  X-36  tailless  aircraft,”  J. 
of  Guidance,  Control,  and  Dynamics,  Vol.  24(5):  903-909, 
2001. 

[3]  A.  J.  Calise,  R.  T.  Rysdyk,  “Nonlinear  Adaptive 
Flight  Control  using  Neural  Networks,”  IEEE  Control  Sys¬ 
tems  Magazine,  pp.  14-25,  Dec.  1998. 

[4]  A.  J.  Calise,  S.  Lee,  and  M.  Sharma,  “Development 
of  a  reconfigurable  flight  control  law  for  the  X-36  tailless 
aircraft,”  J.  of  Guidance,  Control,  and  Dynamics,  Vol.  24(5): 
896-902,  2001. 

[5]  J.  A.  Farrell, “Stability  and  Approximator  Conver¬ 
gence  in  Nonparametric  Nonlinear  Adaptive  Control,”  IEEE 
Trans,  on  Neural  Networks,  Vol.  9(5):  1008-1020,  1998. 

[6]  J.  A.  Farrell,  M.  Sharma,  and  M.  Polycarpou,  “Longi¬ 
tudinal  Flight-path  Control  using  on-line  Function  Approx¬ 
imation,”  AIAA  J.  of  Guidance,  Control  and  Dynamics,  ac¬ 
cepted  for  publication  Jan.  2003. 

[7]  E.  Johnson  and  A.  Calise.  “Neural  network  adaptive 
control  of  systems  with  input  saturation,”  In  American  Con¬ 
trol  Conference,  June  2001. 

[8]  M.  Krstic,  I.  Kanellakopoulos,  and  P.  Kokotovic,  Non¬ 
linear  and  Adaptive  Control  Design,  New  York,  John  Wiley, 

1995. 

[9]  S.  P.  Karason  and  A.  M.  Annaswamy.  “Adaptive  con¬ 
trol  in  the  presence  of  input  constraints,”  IEEE  Trans,  on 
Automatic  Control,  Vol.  39  (11):  2325-2330,  1994. 

[10]  R.  Monopoli.  “Adaptive  control  for  systems  with  hard 
saturation,”  In  IEEE  Conf.  on  Decision  and  Control,  pp. 
841-843,  1975. 

[11]  M.  Polycarpou  and  P.  Ioannou,  “A  Robust  Adaptive 
Nonlinear  Control  Design,”  Automatica,  vol.  32(3):  423-427, 

1996. 

[12]  M.  Polycarpou  and  M.  Mears,  “Stable  Adap¬ 
tive  Tracking  of  Uncertain  Systems  Using  Nonlinearly 


(a)  Initial  state  and  control 


(b)  Final  state  and  control 


Figure  1:  Plots  of  the  aircraft  state  and  control  signal 
at  the  beginning  and  end  of  a  150  s  simulation 
with  7 r  £  ±5  deg:  command  (dotted),  filtered 
command  (solid),  actual  (dashed). 

Parametrized  On-Line  Approximators,”  Int.  .1.  of  Control, 
vol.  70(3):  363-384,  1998. 

[13]  M.  Polycarpou,  J.  A.  Farrell,  and  M.  Sharma,  “On- 
Line  Approximation  Control  of  Uncertain  Nonlinear  Sys¬ 
tems:  Issues  with  Control  Input  Saturation”  In  Proc.  of  the 
2003  American  Controls  Conference. 

[14]  R.  T.  Rysdyk,  Adaptive  Nonlinear  Flight  Control, 
Ph.D.  Thesis,  Georgia  Institute  of  Technology,  School  of 
Aero.  Eng.,  Atlanta,  GA,  1998. 

[15]  S.  N.  Singh  and  M.  Steinberg,  “Adaptive  Control  of 
Feedback  Linearizable  Nonlinear  Systems  with  Application 
to  Flight  Control,”  .1.  of  Guid,  Control,  and  Dynamics, 
19(4):  871-877,  1996. 

[16]  B.  L.  Stevens  and  F.  L.  Lewis,  “Aircraft  Control  and 
Simulation,”  John  Wiley,  1992. 

[17]  H.  Wang  and  J.  Sun,  “Modified  model  reference  adap¬ 
tive  control  with  saturated  inputs,”  In  Conf.  on  Decision  and 
Control,  pp.  3255-3256,  1992. 

[18]  D.  G.  Ward,  J.  F.  Monaco  and  M.  Bodson,  “Devel¬ 
opment  and  flight  testing  of  a  parameter  identification  al¬ 
gorithm  for  reconfigurable  control,”  AIAA  J.  of  Guidance, 
Control,  and  Dynamics,  Vol.  21(6):  948-956,  1998. 


Table  1:  Simulation  Parameters:  The  first  row  of  cells  defines  the  parameters  of  the  model.  The  second  and  third  rows  of 
cells  define  controller  parameters.  The  fourth  row  of  cells  defines  the  state  constraints  and  saturation  limits. 


Parameter: 
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Let 

M0 

Mq 

Ms 

Value: 

-0.1 
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r5 

Value: 

0.4 

16.0 
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30.0 

Parameter: 

k1 

kQ 

Value: 

1.3 
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(a)  Tracking  errors  x. 


(b)  Tracking  errors  x. 


Figure  2:  Plots  of  the  aircraft  tracking  errors  for  7r  £  ±5  deg. 


(a)  Lyapunov  quantities. 


(b)  Parameter  estimates. 


Figure  3:  For  7,.  £  ±5  deg:  Left  top  -  Plot  of  y  0(t)T T_10(f)  versus  t.  Left  middle  -  Plot  of  ||a:||  versus  t.  Left  bottom  - 
Plot  of  a yV(t)  versus  t.  Right  -  Plots  of  each  component  of  0  versus  t. 


Lyapunov  norm  Norm  of  state  error^  Norm  of  parameter  error  V'de9/s  ^  <Wde9  ^  Wde9  E  ,  ,  ,  Q’de9/S  ,  ,  Y’de9 


e  4:  Plots  of  the  aircraft  state  and  control  at  the  beginning  and  end  of  a  150  s  simulation  for  £  ±10  deg:  magnitude 


limited  command  (dotted),  filtered  command  (solid),  actual  (dashed) 


Figure  5:  Plots  of  the  aircraft  tracking  errors  for  7 r  £  ±10  deg. 


(a)  Lyapunov  quantities.  (b)  Parameter  estimates. 


Figure  6:  For  £  ±5  deg:  Left  top  -  Plot  of  y  0(t)TP_10(t)  versus  t.  Left  middle  -  Plot  of  ||a:||  versus  t.  Left  bottom 
Plot  of  a JV(t)  versus  t.  Right  -  Plots  of  each  component  of  0  versus  t. 


